Periodic orbit theory and spectral statistics for scaling quantum 
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The explicit solution to the spectral problem of quantum graphs found recently in 
2(|, is used to produce the exact periodic orbit theory description for the probability 
distributions of spectral statistics, including the distribution for the nearest neighbor 
separations, s n = k n — k n -\, and the distribution of the spectral oscillations around 
the average, 5k n = k n — k n . 
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I. INTRODUCTION 

Quantum graphs consist of a quantum particle moving on a quasi one-dimensional net- 
work. In the limit h = 0, quantum networks produce a nonintegrable classical counterpart - 
a classical particle moving on the same network, scattering randomly on its vertexes jl[ • 
As shown in this stochastic dynamics is characterized by an exponential proliferation of 
periodic orbits, positive Kolmogorov entropy and other familiar features of finite dimensional 
deterministic chaotic systems. 

Such classical nonintegrability is clearly manifested in quantum regime. Quantum net- 
works provide excellent illustrations to many general concepts, phenomenological hypothe- 



ses and mathematical constructions o 



studies fly 



and analytical 



3H,y,Q. 



quantum chao 
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ogy. For example, extensive numerical 
15J1 have demonstrated that the statis- 



tics of the nearest neighbor spacings distribution, the two point autocorrelation function, 
the form factor, and the spectral rigidity of the quantum graph spectra are close to the 
ones predicted by the random matrix theory (RMT). Since the latter three statistics can be 
expressed in terms of the spectral density functional, they were also studied analytically in 



terms of the Gutzwiller's periodic orbit series expansion , Q, Q, Q, Q, Q, H, Q . 

Due to the relatively simple structure of the network dynamics, the results of the periodic 
orbit theory analysis of quantum graphs are particularly complete. Moreover, the periodic 
orbit expansions which usually have semiclassical accuracy, are exact for the quantum net- 
works and can be viewed as mathematical theorems. In addition, as it was shown recently 
m , the periodic orbit theory for the quantum graphs can describe not only the 

global characteristics of the spectrum (e.g. the density of states, spectral staircase, quan- 
tum and classical dynamical zeta functions, etc.), but also the individual eigenvalues of the 
energy or the momentum. 

This fact provides an interesting opportunity to study several additional statistical dis- 
tributions, including those that are not directly accessible via the Gutzwiller's expansion for 
the density of states, such the distribution of the eigenvalue fluctuations around the average, 
b n = k n — k n , the nearest neighbor spacings, s n = k n — k n -i, etc., which is the main subject 
of this paper. 



The paper is organized as following: Section II reviews the spectral hierarchy method 20]. 
Section III discusses the statistical spectral distributions for the regular quantum graphs, 
which are later generalized for the irregular graphs in Sections IV and V. A short discussion 
of certain statistical universality aspects of the resulting distributions is given in the Section 
VI. 

II. SPECTRAL HIERARCHY FOR QUANTUM NETWORKS 

The idea of producing the individual momentum eigenvalues k n = kn^ is based on using 
the periodic orbit expansion for the density of states, 



p (0) 



oo 

(A0 = £<^-£i o) ), (i) 



n=l 



and an auxiliary sequence kn , that separates the spectral points from one another: 

k n 1 l 1 <k£ ) <k n 1) , n = l,... . (2) 
From these two constituents one can obtain the quantum energy levels via 

fc£°>= / p<®{k)kdk. (3) 



3 



If any sequence with property ([2]) is known as a global function of n, kn^ = fcW(n), the 
relationship (J2J) produces the explicit solution to the spectral problem, 

The explicit integration in ([3]) is possible due to the exact periodic orbit expansion for the 



density of states. As shown in 



the exact expansion for p^(k), has the form 



01 EL 

p (o) {k) = h + I Re VLfAfe^ ^. (5) 

P 

where Lp " 1 and are correspondingly the action length and the weight factor of the 
periodic orbit p, L is the total action length of the network. For the scaling quantum 
graphs the weight coefficients, Ap , are independent. 

Effectively, obtaining the spectral points as a function of their index is equivalent to 
"inverting" the spectral staircase function, 

N(k) = Y^e(k-k i ) 1 (6) 

i 

i.e. passing from N(k n ) = n to k(n) = A^ _1 (n). Geometrically, finding an auxiliary sequence 
that singles out separate peaks in p]), amounts to finding a suitable monotone function f(k), 
whose graph intersects every stair step of the spectral staircase. As it is illustrated on FigJH 
the intersection points, kn\ 

f (W) = N = n, (7) 



clearly satisfy the condition (T5]), so solving the equation (J7J) would yield the whole sequence 
kn as a single globally defined function of the index n. It is well known however, that 
following the behavior of the spectral staircase function ([6]) in such detail is generally a 



difficult task (see e.g. |23|). Luckily, quantum graphs allow an alternative approach, [17, 
I20I ]. which is based on certain properties of their spectral determinant. 
As shown in (2), Q, Q, [2^] , the spectral determinant A(/c) for quantum graphs is a finite 



order exponential sum, 



A(fc) = 1 + e* 2 ^-^) - a ie i2(Llfc -" 7l) , (8) 
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FIG. 1: A continuous function f{k) that pierces every stair step of the spectral staircase N{k) of 

an 4 vertex linear chain graph. The intersection points kn^ separate the momentum eigenvalues 

L(o) 

with constant cij, Lj < L and 7$. The order iVr of the sum depends on the topology of the 
graph T. The explicit form (jSj) can obtained from imposing the boundary conditions on the 



wave function o 
quantization 0, 



the quantum particle moving on the networ 



£, e.g. by using the scattering 



1?| or the Bogomolny's transfer operator 24| methods. The roots of A(k) 



define the quantum spectrum of the momentum, A(kn^) = 0. 

There are three key properties of the spectral determinant A(k) relevant 
discussion. First, its roots as well as the roots of all of its derivatives are real 



or the following 



31 



321 ] . Second, 



there is exactly one root of its jth derivative, A^(k), between every two neighboring roots 

Am 
32J. 



31 



of A" +1 '(A;). This implies that the zeroes of A(k) are interlaced by the zeroes o: 
as required by (J2J), which in turn are interlaced by the zeroes of A"(k) and so on 
Last,, a, shown in Q, the hi g he r is the order of the derivative, the more orderly is the 
behavior of the roots of A^\kn ) = (Fig. [2]). In fact, for any quantum graph system, 
there exists a finite integer r (called the regularity degree of the graph in [2^], see also 25],) 
such that the roots of A^ r \k) can be interlaced by a periodic sequence of points 



According to the criterion used in 
for which the inequality 



17 



18 



4) 



(9) 



20j ] , the regularity degree r is the minimal integer 



E 



Li 
L 



< 1 



(10) 
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FIG. 2: Histogram of the fluctuations of the roots of the 6-star graph of irregularity degree 5 and 
of the following 5 separating sequences, S n ° , Sn\ The initial spread of the fluctuations 8n^ 
becomes progressively narrow for the higher values of j. 



holds. Since (j^j < 1; this criterion allows to find a finite regularity degree for every 
quantum graph system using the coefficients of the spectral determinant ([8]). Hence, there 
exist r + 1 almost periodic sequences of interest, kn\ j = 0, r, such that 

A (i) (ty) = 0, (11) 



(12) 



and 

The density functional for the jth sequence, 

n 

allows to pinpoint the exact location of kn on the interval between kn and &i J 2i, 



(13) 



K n-1 



p {j - 1] {k) kdk, 



(14) 



for j — r + l,r — 1, 1 and n = 1,... . Since the roots of each A^(/c) form an almost 
periodic set [32| the density functional of each of the sequences kn can be expanded into 
an explicit harmonic series, which allows to evaluate the integral f[T4|) explicitly for every n 
and j. 

The strategy that allows obtaining the separating sequences kn and eventually getting 
the physical spectral sequence kn \ follows directly from the bootstrapping (i.e. interlacing) 
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FIG. 3: The bootstrapping of the spectral staircase (dotted line) by the N^ 1 ' staircase (dash- 
dotted line), bootstrapped in turn by the staircase (solid line) intersected by the Weyl's 
average (straight line). Note that the staircase does not bootstrap the N^°\ These graphs 



were obtained for the dressed 4 vertex linear chain graph 



1 



t(r) 



property of the sequences ( fl2l) and the relationship (Il4"l) . One first finds the sequence k 
starting from the periodic separators ©, then uses it to find kn , and so on. After r steps 



of moving up the hierarchy the spectrum kn^ is produced [20|. 



Geometrically, this algorithm can illustrated using the integrated densities 

N®(k) = ( k ~ ~ k n ] ) • ( 15 ) 

n 

The separating property of the jih sequence with respect to the sequence fcn implies 
that the staircase N^'(k) "bootstraps" the staircase N^(k), whereas N^(k) bootstraps 
N^(k) and so on (see Figj3l compare to Fig{T]). The final staircase N^(k) is pierced by 
the Weyl's average 

N{k) = ^k-\. (16) 

7T Z 

In the simplest case of r = 0, the spectral points themselves can be separated from one 
another by a periodic sequence (Q. Geometrically, this implies that the Weyl's average 
pierces every stair step of the original N(k), which guarantees the existence of periodic 



separators (Q. Such networks were referred to as the "regular graphs" in 



17 



18 



2fJ. The 



result of integration ([3]) in this case produces the quantum eigenvalues in the form of the 



7 



periodic orbit series, 

k n = ^-n - — -m sin ( -y- sin K° )n ) ' ( 17 ) 



where uf^ = nLp /Lq. The first term in f|T7|) gives the average behavior of the eigenvalue 
sequence, k n = j^, and the subsequent periodic orbit sum describes zero- mean fluctuations 
of k n s around the average. It will be more convenient to describe the fluctuations in terms 
of the quantity Sn^ = Lq (k — k n ) /tt, n = 1, 2, 

o AO) /,,(°)\ 
5 n ] = ~ E ^ 8111 \^tj Sin M° )n ) • 



As demonstrated in 17, m 2oi, if the sum (JT7|) includes only the orbits that involve a 
certain fixed number m of scatterings at the vertexes of the graph, it produces the mth 
order approximation to the exact value of k n for each n. 



III. EIGENVALUE DISTRIBUTION FOR REGULAR GRAPHS 

Let us first study the spectral fluctuation statistics for the regular graphs based on the 
expansion (ITT)) . The transition to statistical description of the sequence 5n \ n = 1,2,..., 
can be made based on the properties of the sequence of the remainders 

X n = [Hmodl> U = i' 2 '-' ( 19 ) 



It is a well known number-theoretic result |33|, [34| that the sequence (Tl9l) is uniformly 
distributed over the interval x n e [0, 1] for any irrational number a. 
For every periodic orbit p in (fl8l) . the frequency oo^ is defined as 



f = m^ih + ... + m% B n NB = (m£\n p ) , (20) 



where fli = tiU/ Lq are the "bond frequencies" defined by the bond lengths k, fij = 1, and 
the vector mj, ' = (m^\, ...,rrip NB ) gives the number of times the orbit p passes over the bond 
i. Considering that for every the function sin(a;p ^n) can be reduced to a combination 
of basic harmonics of f^n, one concludes that in the generic case in which every Zj / Lq is an 
irrational number, the phases 

x i,n = [^n] mod2n . (21) 
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in each term in (|20|) will generate random outputs, uniformly distributed in the interval 
[0, 2n}. Hence, in the context of studying the statistical properties of the eigenvalue sequence, 
the argument of every factor sm(uj^n) in ( TTTj) can be treated as a function 

sin (uj^n) — > sin (mp\xi + ... + m p % B XN B j = sin (rh^x) , (22) 

of Nb — 1 independent random variables Xi, that are distributed in the interval [0, 2n]. Hence 
the set of the deviations of the eigenvalues from the average is statistically described by a 
series or random inputs corresponding to the periodic orbit expansion f|T8|) 

o 4(0) 

^ 0) = ~ E sin \~t ) sin • (23) 

The maximal amplitude of an input corresponding to a periodic orbit p coincides with the 
amplitude of the same orbit's contribution into the exact periodic orbit sum (p3 
It is now a straightforward task to obtain the distribution of Sn^ via 

N B -1 



pf (*«») = jT ... jT * ^ (0) + E C T sin K 0) -) J II ^- ( 24 ) 



where 

2 a<°» / w (0 > 



Using the exponential representation of the ^-functional, one has 



P f (tfW) = / dke^ff (jfe, cf ) , (26) 



where the characteristic function (k, Cp\ is defined explicitly via the periodic orbits 
and graph parameters, 



Nh—1 



(jfe,Cj°)) =/"..•/ V^^K'*) "Q (27) 



As in the case of the series expansion for a finite order (rath) correction to the exact result 

(0) (0) 



is obtained by considering the orbits that involve the same number \rrip \ = ra^ ] + ... +rn^' 



of vertex scatterings. 

In general, the statistical properties of a spectral quantity z n that has a periodic orbit 
series expansion 

^> = /f -E 4° W4 0) * ( 2 8) 
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where ip^ is an n-independent phase and ff^ 1 is a shift term (see below), are described by 
the random sequence 



4 0) 



/ 2 (0) -E4 0)cos (4 0) ^+^ 0) )- ( 29 ) 



The corresponding characteristic function of probability distribution for z will have the form 



2?T /-27T 



Fj®(k) = I ... I e iki: ^ 0)coB ^ 0)s+ ^ 0) )dx, (30) 







where dx = Ui=i f 1 - 

Another case of spectral characteristics that can be treated by this approach is provided 
e.g. by the separation between two eigenvalues, k n+m and k n , for a fixed m. This quantity 
is defined by the expansion 

<C = f - " E D fk COS ( W W„ + ^) , (31) 



with 



4 Af> /ui 0) \ (J®m 



Ba=^sm(^|sin|^], (32) 
which gives rise to the random series 



<e = ™ - e ^ o) c ° s ( 4 o) - - ) . (33) 



V 



The distribution of sS' can be obtained as 

^00 = jdkA^FjV^Dfl) (34) 

where Fj$ (k, D$ty is defined according to ( l30il . In the case m = 1, these formulae de- 
scribe the nearest neighbors distribution. Clearly, for the difference between the fluctuations 
themselves, <Jn°m = 5 n+m - 5 n , the distribution is P$l{o) = P$ + ^J. 

One could also study the mean of the two neighboring deviations, ffl = + d^lA /2. 
This characteristics can be used to generate the separators directly from kn^ , kn = n + ffi 
independently of the properties of the spectral determinant. The corresponding series has 
the expansion coefficients 

4»» = i^|si„K«)) (35) 
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(0) u/ 0) 

and phases (p p = | 1-. The energy fluctuations, 5E n , have the expansion coefficients 

H > nw[^ sm [—)- ms — 1 (36) 

phases (pf^ 1 = and f E = etc. 

These quantities can be used to find the periodic orbit expansions for higher order statis- 
tics, such as the correlator /^n^i+ny or the autocorrelation function, 

N 



R2{X) = %{p{ k + l)p{ k -l)) = Tj^T,Y, 5 -K + x). (37) 



-2 

J ° ' ~ ~" n=l m^O 

which defines the probability to find a new level at a distance x ^ from a given old one, 
whether or not these levels are nearest neighbors. The Fourier image of R,2(x), the form 
factor, K 2 (r) is given by 

N N I \ 

K 2 (t) = — lim - VV e i(kn +m -k n )r = JL(y^ e -i{ku +m -k n )r \ _ , 3g N 

n=l m=l \ m ri 



For the regular graphs with k n given by ffTTj) . the averaging over n can be performed via 

N 

k 2 (t) = — iim lyye^-^*™^^!)))^ 

n=l m 

= —J2 [ e ! ^^ cos K , ^°T) Yl —■ (39) 

The latter integral is defined by (l3lj) . so 

oo 

^(r) = fE e " !fTf S(^S' (40) 



m=l 



Formula (HOI can also be obtained directly from (1551) by averaging over the random variable 
Sm using the distribution (IMj) . 



\ m l n m 

Hence R2{x) is given by 

^(*) = ££/ ^-^S M20 = t>(*). ( 41 ) 
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as the sum of the probabilities that the two eigenvalues separated by the interval x have 
m — 1 other eigenvalues in-between. 

It should be emphasized that all the probability distributions above are obtained in 
the context of the standard periodic orbit theory framework. All the distributions for the 
regular level fluctuations derived in this section are closed, self contained expressions, defined 
in terms of the periodic orbits and graph parameters. It is also important to notice that 
the statistical properties of some (especially some regular) quantum graphs, despite being 
strongly stochastic in the classical regime, deviate from the universal Wignerian distributions 
predicted by the RMT. However these well as the irregular graphs discussed below, 

are equally well described via statistical description of the periodic orbit expansion series 
for the spectral sequences. 



IV. SPECTRAL EXPANSIONS FOR IRREGULAR GRAPHS 

As mentioned in the Section II, one can find the roots of /S.^~ l \k) by using the density 
p^\k) and the separators fc„ in formula (fT4j) . which yields 

fc(i-D = kU) N u-D _ k^ iN u-D _ r N^\k)dk. (42) 



n-1 



Every staircase function N^\k) can be decomposed into the average and the oscillating 
parts, N®(k) = NW(k) + SN^(k), where the average integrated density for every j is 

N^(k) = —k- - = N(k). (43) 

7T 2 

The bootstrapping flT2} of Ui~ X) by k { n j) (or N^~^(k) by N®(k), see Fig|3J implies that 

N (j ~ 1} = n. (44) 

Using (HHD, ( 1441 and writing kn^ in the form 

fcW = 1- (n + 6$) , (45) 
we get for the j — 1 generation of the separators 

f Ci) 

kt 1] = t + It (* -) - - lr ( {&2 - ( *- )2 ) - IJ 5iv0 '" 1)(A;) (46) 
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The harmonic series expansion for 5N^(k), 

SN^(k) = -lmJ2 Afe lL ? )k , (47) 
v 

allows to compute the integral in f|4*8|) explicitly, which yields the fluctuating part of kn 
via an expansion similar to (Fl8|) . 

st i] = ft l] (ee*) -Y, c t i] (ee) - (-rw-) (ee)) ,m 



where now the expansion coefficients, 



2 A 



ct i] = sin 4- (*» - + x ) ' (49) 

Up 

the "zero shift" term, 

ft* = \ (#> - <&) - \ (OW - (€!) 2 ) , (50) 



and the phases, 

r(i) _j_ jtC?) _ i 

< j = o < \ (51) 



are functions of the fluctuations 5„ and b^-i on the previous level of the hierarchy. In the 
particular case when r = 0, 8n = = 1/2, ( HHl) coincides with the oscillating part of 

The equation (14"5|) shows that the fluctuations 5n ~ have two sources. In addition to the 
oscillations induced by the periodic orbit sum in (jUJ), there are also oscillating contributions 
produced by 5n and that bring in the oscillations from all the previous levels of the 
hierarchy. This equation will later be used to produce the exact statistical distribution for 
the fluctuations 8$ at each level of the hierarchy. 

As an example of where the coefficients A® can be o 

a simple example of a 2-bond regular graph, discussed in 
graph is a strictly regular system, which does not require auxiliary separating sequences 
for obtaining its spectrum, it is nevertheless useful to use this case to have an immediate 
illustration of the explicit form of the coefficients Ap\ 



- 1 "btained directly, one can use 
Q. Although the 2-bond 



As shown in 



17 



IS 



21| , the spectral equation in this case has the form 

sm((h + l 2 )k) - r sin((Z 1 - l 2 )k) = 0, (52) 



13 



where r is the reflection coefficient at the midd 
engths. This equation was obtained in 



18 



e vertex and l\ and I2 are the two bond 



21| using scattering quantization method 



l2|] as A(fc) = det(l — S(k)) = 0, where scattering matrix S(k) = S^(k) in this case is 



S®(k) 



( 







\ 










te il 2 k 


te ihk 










{ 







) 



(53) 



The unitarity of S(k) is guaranteed by the "flux conservation" relationship between the 
reflection and the transmission coefficient, t 2 + r 2 = 1. The expansion of the spectral 
determinant , yields the exact periodic orbit expansion for N(k), 



N(k) = N(k) - -hn^2 ^ [(-1)^)^)^]" e iuLpk 



(54) 



p.y 



where v is the multiple traversal index for the prime periodic trajectory p, a(p) and 2r(p) 
are the numbers of reflections and transmissions for p at the middle vertex, and the factor 
(— 1)x(p) defines the Maslov index 21|. The action length of p is L p = m p> ili + mp^h, 



according to the number of times, m p 1 and m P 2, it traverses the bonds l\ and ^2- One can 



notice that the equation for the separating points k 



(i) 

n , 



cos((/i + l2)k) — rcjcos((Zi — h)k) = 0, 



(55) 



where u — (1% — h)/(h + ^2) < 1, reduces to the original equation (|52"]) . if r — > = cjr and 
Zi/c — > Zi/c + ~. Therefore, the expansion for N^(k) can be obtained from det(l — ^■'(/c)), 
where unitary matrix S^(k) is obtained via a two parameter deformation of the unitary 
matrix S^°\k), 



( 

r (i) e «i*+f 



— e* 



V 



£(l) e «l 


















j-(l) e ihk 
_ r (i) e ihk 





\ 



(56) 



with = tor and t« = V 1 - (r«) 2 . The weight coefficients are now explicitly 
defined via the products of the matrix elements of S^(k), just as the coefficients were 
defined via S^(k) in (El). 
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Clearly, all even order derivatives of the spectral equation (152]) for j = 2w have the 
form sin((Zi + l2)k) — rto 2w sin((Z 1 — l2)k) = 0, so with the replacement r^ 2 ^ = cu^^r, 



t( 2w ) — y 1 — (H 2 "")) 2 , the form of the coefficients A p 2w ^ is structurally same as the one 
found in ( 1541) . The odd degree derivatives A p 2w+1 ^ are produced by the expansion of the 
determinant det(l — S^ 2w+1 \k)) analogous to ( 156]) . It should be mentioned however, that in 
general the task of obtaining the exact form of the coefficients A p ^ and frequencies ujp^ for 
j 7^ is not as straightforward as in this simple case and requires a more detailed analysis. 
In particular, the harmonic expansions such as (148]) may include bond combinations that do 
not correspond to connected periodic orbits. 

Using the expansion for the fluctuations Sn ^ , one can find the harmonic expansions for 



other spectral characteristics. For example, the m-neighbor difference, Sn,m^ = m + S^ +n { 
Cj-i) 



Sn , the expansion is 

,Cj-i) _ f(j-i) (JJ) Jj) JJ) rfjWVD^'ros fnw^ + J^ 1 ') (^7) 



where 



f(j-i) _ Jj) 1 I (Jj) _ Jj) \ (om - J j) - s (i) ^ 4- /V j) - ^ (hti) 
Js s n,m* 2 \jn,m s n,m-lJ °n,m s n-l,mJ T Sn \^n— l,m °n,mj \ ua ) 

and Cn 0) = I + 5^2iJ- The harmonic expansion coefficients in this case are 



4(J-1) , , ,(3-1) , ,0-1) 

- 12 ( S i n 2 ^ JJ) + sin 2 ^ 

Ulp - - 



,(j'-i) , ,Cj-i) \ , 



2 ^Sin ^ S^ m Sin ^ S n-l,m^ COS ^ ^ S njm-1 S r^m S n-l,m^ ) (59) 

In case m = 1 the system of equations (j57p yields the nearest neighbor distances between of 
the j — 1 level separators. 

The nearest neighbor average, $\ that appears e.g. in (159]) has the expansion 

where 

_ 1 L| / (i) (i) _ 2tt\ / U) _ 2tt 2tt (j) 
; « " 4 7r 2 V S " L ) V ™ L A) 1 

r(3-i) _ sin=^ (s {j) + s {j) 
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FIG. 4: Histogram of the nearest neighbor separations Sn^ (bottom distribution), Sn and 
(top distribution) for the fully connected quadrangle graph (top right corner) of the irregularity 
degree 3, obtained for 80,000 roots of the corresponding spectral equations. The maximal nearest 
neighbor separation for j = in this case is s max = 8.68, while regular cell size is tt/Sq = 2.28. 

The expansions of the kind (JHJ), floTj) and (160]) . provide an exact description of the prop- 
agation of the spectral characteristics across the hierarchy in terms of the geometry of the 
graph. For the graphs of high degree of irregularity, the equations ( )48l) . ( 1571) . ( |60j) . etc., can 
be considered as discretizations of nonlinear differential equations for continuous functions 

~ S(n,j), sffm ~ s m (n,j), ~ £(n,j), etc. 

V. FLUCTUATION STATISTICS 

As demonstrated in the previous section, the fluctuations Zn (e.g. (1481) . fl57|) . (160]) ) at 
the jth level of the hierarchy depend on the fluctuations on all the previous levels as well as 
on the oscillations introduced by the harmonic terms at the level j, 

= f«-*> (4 i} ) - E ci r l) «>) cos ("t 1)n + , (6i) 

As in the regular case considered in Chapter III, probabilistic description of the sequence 
z^ at each j is obtained by considering as a function of a random vector x, generated 




16 



by the sequence x { = [^n] mod27r , 

zt 1] = /°'- X) (*?>) - E c ?" } (4°) «* (rh^x+^-V (,0))) . (62) 

p 

Since the index n and the "bond frequencies" Qi are the same for all the z£* expansions 
across the hierarchy, the random variable x is the same in all the harmonic expansion terms 
in the stochastic series fl62l) . j = 0, r. This implies however, that unlike the regular 
case, the harmonic oscillations at the higher levels of the hierarchy (j < r) are not the only 
source of randomness in the series ( 1621) . Additional "noise" at the level j is injected into 
the equation ( 1621) by the variables z x , that introduce the fluctuations from all the previous 
levels of the hierarchy into its j — 1 level. The expansions ( 1621) allow to obtain the probability 
distributions for each of the z^'s from 

P 0--i)( z a-i)) = J § ^(i-i) _ Z U-D (4.) ^CH-D (...), x ) , x )) dx , (63) 
where dx — Y[ ip- The "nested" structure of the Zx in ( 1631) can be unfolded in the form 

i 

pO'-D^O-D) = j 5(z^-z^(z ( i\x)) 

•■•^(^-^))nn^- ^ 

j « 

in which the chain of the delta functionals above can be understood as the conditional 
probability densities for obtaining the value z^~ x \ given the values of x and 
z^ . Here the index i runs over the number of the elements zf^ that appear in the expansion 
fl6Tj) of the corresponding quantity z^~ x \ For example, the expansion of 5n depends on 
5^2 1 and Sn \ so in this case index i takes two values, i — 1,2. In the following this index 
will be omitted for conciseness of notations. 

The hierarchical organization of the fluctuations suggests a natural way of approximate 
evaluation of the expression (1641) . Since for every j the function z£* is a complex, rapidly 
oscillating function of x, and since every z^" 1 ' characteristics depends explicitly only on 
the previous level sequence z^\ it is natural to consider a simple approximation to (I64"I) . 

in which z^> in the equations (l62l) are treated as independent random variables distributed 

(i) 

according to P 2 , 

P 0'-i)(^0'-i)) = /" tf ^(j-i) _ z U-i) ( z <j) iX )) P^dz^dx. (65) 
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A more formal and elaborate argument that will not be discussed here in detail is based 
on inverting the dependence z^~ x ' = z^~ l \x) for each j in the arguments of the delta 
functionals in (164j) (i.e. in effect using Bayes relationships P(z\x)P(x) = P(x\z)P(z)). It 
is clear that due to the complex oscillatory behavior of the expansions z^~ l \x), one can 
approximate each of the P^\x\z) distributions by a uniform distribution, which effectively 
corresponds to introducing separate x- variables for each level j, and leaves, after intermediate 
integrations, the distribution P^ in fl65|) . 

Such approach takes a natural advantage the hierarchical organization of the fluctuations 
produced by the sequences . Since the distribution Pz^ at the regular level can be 
computed directly, the distribution Pj r 1 \ pj r 2 \ can be computed in sequence, from 
the previous levels of the hierarchy to the next. The final distributions P^ will apply to 
the physical characteristics of the spectrum. 

As an example, one can consider the transition from the level j to the level j — 1 of the 

f i) (i) 

distributions for the deviations from the average. The jth level fluctuations Sn and 5^_ 1 in 
( HSjl . can now be considered as random variables, 5± and 62 , distributed according to Pg (5). 
Hence for the density P^ (8) one can write 

Pt l) (*) = / dhdfodxPP (5,) P^ (5 2 )5(5 - ft 1] (8?, 5?) 

+ E C ?~ 1) (5?,8?)sm(m^x + <p^ (5?, 8?)) ) . (66) 
p ' 

Exponentiation of the delta function produces 

P^~ x \8) = J dke ikS (k, C¥~ 1] )) nU) , (67) 

where the function 

Ft 1] (k) = r ... r ^^-^k^-w- 15 )^ m 

Jo Jo 

produced by the integral over XiS is reminiscent of the (1271) . The parenthesis (*, *) n o) denote 
averaging of the F^\k) with the weight 

flW^, j®, k) = e'^^^pU) (5l) P W {8%) ( (69) 

which corresponds to averaging over the "separator disorder" , which yields the characteristic 
function of the distribution (167|) . In case if Pg = 5 [8^ r+l ^ — ~) (ordered separators at 
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j = r), one recovers the regular level distribution (1261) . For j < r, obtaining the probability 
distributions Pg 1 \8) for requires an additional averaging over the disorder produced by 
the fluctuating sequences of the separators kn- 

The expression for other spectral characteristics have similar structure. In the case of the 
m-th neighbor separation statistics, the probability distribution for sli is given by 

Pt 1 H*) = j e ik{s - m) (F^)^dk, (70) 

where the "characteristic exponential" F^} generated by the harmonic series 5sx,rn^ similarly 
to (|68|) . is averaged with the weight 

(*i,* a ,«s,e) = e-t^wtipU \t)P®(8{)P!£(8 2 )pV>(8 3 ). (71) 

In the case if m = 1 this yields the distributions for the nearest neighbor spacings. The 
resulting distributions can be used to compute the form factor, fQ ( r )> the autocorrelation 
pjf\x) function and so on. 



VI. TRIGONOMETRIC SUMS AND SPECTRAL UNIVERSALITY 



An important advantage of obtaining the statistical description of spectra in terms of the 
harmonic expansions, is that the well known universality features of the quantum chaotic 
systems (e.g. 



29 



shown e.g. in 



30| ) can now be analyzed within the context of the theory of trigonometric 

35 . 3(1 and references therein). As 



series and weakly dependent random variab 



37 



40 



41 



42 



43. 



44 



45 



es see 



461 ]. under certain conditions separate terms of 



trigonometric (or more general 39]) series statistically behave as a set of weakly dependent 
random variables 53] . This allows establishing analytically certain universal features for the 
distribution of their sums, e.g. the corresponding generalization of the central limit theorem. 



It is well known 
series 



37|, 



41 



42 



43|, 



44 



45 



46] ) , that if the frequencies of the trigonometric 



fi x ) = ^2°k cos(2irn k x + <f> k ) 



increase sufficiently rapidly so that the "critical condition" , 



1 + 



oc 



(72) 



(73) 
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is satisfied, then a central limit theorem holds for the sum (I72p . 

1 - 

— ^ °k cos(2im k x + 4> k ) -> Ao,i, 



A" 



(74) 



k=l 



where A/o,i is the standard normal distribution and = \ ^^ fc=1 c^. 

Although the series (I72p is much simpler than the spectral expansions studied in the 
context of the periodic orbit theory, similar result was previously observed in the physical 



literature. Based on the extensive numerical simulations, it was hypothesized in [30J that 
the fluctuations of the spectral staircase SN(k) of ageneric quantum chaotic system are 
normally distributed with the standard deviation 



9 5Z A P' 



(75) 



introduced in 



30t ] via spectral rigidity, where A p 's are the expansion coefficients (j4"7|) for 



SN(k). In particular, this applies to the quantum graphs, for which there exists an exact 
periodic orbit expansion (j4"7j) for 5N(k). Although the fluctuations of the quantum graph 
spectrum are finite |20j, their distributions are typically well approximated by the Gaussian 
distribution (see FigJ2]and Fig J5] below), just as the finite range P s can approximated by the 
Wignerian distribution |2[] (see FigHJ). 

The explicit expansions (|29|) and fl62l) allow to extend the scope of the "central limit 
theorem for spectral fluctuations" of [30j hypothesis to a much wider range of spectral 
characteristics. For example, the probability distributions of the nearest neighbor average, 



£n , shown in Fig. [5J which can be obtained from (1611 as 



(76) 



with the weight 



(77) 



are also well approximated by Gaussian distribution with the corresponding variance. 
It is important to mention that alt hough t he conditio n (1751) is required for the general 

461 ]. in some cases it is pos sible to 



37 



40 



41 



42 



43 



proof of the central limit theorem (1741) , 
establish the existence of limit distributions if (1731 is violated (see e.g. [45, 



47 



48 



60|). This 



is essential for the periodic orbit expansions, because lacunarity condition (1731 generally may 
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1.8 




FIG. 5: Histogram of the nearest neighbor average for the even levels of the hierarchy of the 
fully connected quadrangle graph of the irregularity degree 7, obtained for 80,000 roots of the 
corresponding spectral equations. The solid line in the background of each histogram represents a 
Gaussian fit with A given by the corresponding sum (|75p . 



not hold for the prime periodic orbit spectrum of the quantum graphs and of more general 
systems. However extensive numerical and empirical evidence points on the existence of the 
corresponding limiting distributions 30 1. 



There are much fewer proved results about the probabilistic behavior o 



5j, 



55 



56 



the multiple 



57 



58 



59 



60) 



trigonometric series such as the spectral expansions (T2"9~l) or (I6"2"l) (see 
and the references therein), specifically regarding the conditions required for convergence 
of their sums to the limiting distributions. Nevertheless, the proposed connection to the 
spectral theory clearly implies the existence of the limiting distributions for such expansions 



30] from the perspective of the 



and points out the origins of the spectral universality 29 
periodic orbit theory. 

The nontrivial role of the transitions between the probability distributions at different 
levels of the hierarchy (such as (!6T]) or (!70|) ) can be seen on the example of the development 
of the nearest neighbor distributions (see FigHJ). Although the distribution at the regular 
level has overall Gaussian shape, the profile of the distributions P s (s) for j > becomes 
progressively closer to the Wignerian form. 
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VII. DISCUSSION 

The method of obtaining the explicit semiclassical expansions for the individual spectral 
points outlined above is based on the possibility to tie the sequence of momentum eigen- 
values, k n = kn , to a certain base, regular sequence kn +1 ^ defined explicitly as a function 
of the index n. As shown in [20], this can be done either directly, as in the case of the 
regular graphs, where the periodic sequence (1751) is used to place the points k n into a system 
of periodic cells, or indirectly, as in the case of the irregular graphs, where a few auxiliary 
sequences, kl%\ are needed to complete the bootstrapping. 

The system of the auxiliary sequences, kn \ that links the base sequence kn +1 ^ and the 
physical spectral sequence kn , together with the rule of transition from the jth to j — 1th 
level, defines the spectral hierarchy. The depth of the hierarchy, i.e. the minimal number of 
the auxiliary sequences necessary to complete the bootstrapping, expresses the complexity 
of the spectral problem with respect to a particular bootstrapping method. 



The spectral hierarchy used in 20j and in this paper is based on using the sequences of 



the roots of the derivatives of the spectral determinant and the base sequence 

= f (n + ~Y (78) 



L V 2 



The index n 6 N that appears explicitly in (1751) is then carried to the 0th level via the 
rule (114)1 . This scheme allows describe the exact evolution of the separating sequences fc>/' 
from the lower to the upper levels of the hierarchy, and to pass on to a general probabilistic 
description of the spectral characteristics, including the ones that not accessible via the 
Gutzwiller's expansion for the density of states. 

The organization of the spectral hierarchy allows one to follow the accumulation of the 
fluctuations from the regular (rth) to the physical (0th) level. Essentially, this method allows 
to unfold the full scale spectral fluctuations in steps, by distributing the disorder across the 
intermediate levels of the hierarchy, and by passing successively from more orderly to more 
disordered sequences. While the starting sequence is perfectly ordered (e.g. (175)) is periodic), 
every other sequence kn \ j < r + 1, is disordered, and the scale of the oscillations increases 
as one passes from the level j to j — 1 [20J . The probability distributions at every level are 
obtained by averaging over the fluctuations introduced by the periodic orbits at that level, 
as well as over the disorder 6^' inherited from the previous level of the hierarchy. 
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The spectral expansions may also provide a physical understanding of the origins of the 
spectral statistical universalities based on the periodic orbit theory. It is well known that 
under certain conditions (e.g. the lacunarity condition (1751) ). separate terms or specially 
combined groups of terms of the trigonometric series behave as weakly dependent random 



variables 



35 



36]). This allows to establish standard universal asymptotic distributions for 



their sums, in particular convergence to a Gaussian distribution with a certain specific vari- 
ance. Interestingly, the same variance was already conjectured for the universal probability 
distribution profile for the 5n^ fluctuations of a generic quantum chaotic systems [3). 

In addition, in the proposed approach, the propagation of the fluctuations through the 
hierarchy and the build-up of the distributions via the corresponding number of aver- 
agings over the disordered sequences hn^ can lead to the appearance of other universal (e.g. 
Wignerian) profiles, as it is the case for the nearest neighbor separation statistics (see Fig. 

H}. 

On the other hand, it is also important that this approach does not overlook the individual 
features of a particular system for the sake of broad universality, and can provide detailed 
description of the distributions that deviate from the universal behavior (as in the case of 
the regular graphs), as well as the degree of such deviation. 

The work was supported in part by the Sloan and Swartz Foundation. 
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